Saturable discrete vector solitons in one-dimensional photonic lattices 
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Localized vectorial modes, with equal frequencies and mutually orthogonal polarizations, are in- 
vestigated both analytically and experimentally in a one-dimensional photonic lattice with saturable 
nonlinearity. It is shown that these modes may span over many lattice elements and that energy 
transfer among the two components is both phase and intensity dependent. The transverse electri- 
cally polarized mode exhibits a single-hump structure and spreads in cascades in saturation, while 
the transverse magnetically polarized mode exhibits splitting into a two-hump structure. Experi- 
mentally such discrete vector solitons are observed in lithium niobate lattices for both coherent and 
mutually incoherent excitations. 



PACS numbers: 42.65.Tg, 42.82.Et, 42.65.Sf, 63.20.Pw 

Solitons or stable strongly localized nonlinear struc- 
tures, which can elastically interact with linear waves 
and other solitons, have been studied in various systems 
in nature, ranging from astrophysics [l| and ocean waves 
[H, down to Josephson junctions [3] and nanowires 
These localized structures exist due to an exact balance 
between two or more counteracting effects such as, for 
example, dispersion and nonlinearity in the temporal do- 
main In the optical domain, solitons may exist in 
specific materials, such as Kerr and photorefractive ones 
[1,01. On the other hand, solitons occur in different forms 
like incoherent, discrete, and vector solitons, which are 
not directly related to a particular material [5]. Vec- 
tor solitons are composite structures that consist of two 
or more components which mutually self-trap in a non- 
linear medium. Importantly, the individual components 
decay in isolation. The existence of vector solitons was 
first suggested by Manakov in 1974 [8[- Later on, vec- 
tor solitons have been, for example, studied in carbon 
disulfide [9] and photorefractive crystals [10|. In peri- 
odic nonlinear systems, the so-called discrete solitons ex- 
ist due to the balance between nonlinearity and discrete- 
ness [11]. They have been observed in diverse physical 
configurations such as biological systems, charge- transfer 
solids, Josephson junctions, micromechanical oscillator 
arrays, and photonic lattices. However, nonlinear optics 
has overtaken a primacy in the nowadays soliton related 
research [12|, LL3| due to a rather mature technology of 
photonic lattice fabrication. Such photonic crystals have 
periodic distributions of the refractive index and light 
propagation is associated with allowed bands and forbid- 
den gaps, analog to propagation of electrons in crystalline 
lattices Q. 

One-dimensional (ID) discrete vector solitons (DVS) 



originating from the first band have been already in- 
vestigated both analytically 15, 16[ and experimentally 



P3] in nonlinear cubic waveguide arrays (WA). The two- 
dimensional case was studied, too [18[. Finally, multi- 
band vector solitons, in which individual components 
stem from different bands, were also recently suggested 
and demonstrated [19[. The aim of the present study is 
to investigate DVS in media with saturable nonlinear- 
ity. Prime examples of such media are photonic lattices 
in photorefractive crystals. Additionally, semiconductors 
at higher light intensities also exhibit saturation [20[. It 
has been shown that saturation, which may occur in a 
cascade manner in discrete systems [2l|, is responsible 
for the existence of multiple zeros of the Peierls-Nabarro 
potential, leading to free steering of large amplitude soli- 
tons, and stable propagation of inter-site modes in ID 
and 2D systems [22|. Various species of two-component 
saturable DVS have been investigated recently, where it 
was assumed that both components have the same po- 
larization and different frequencies [23]. In what follows 
we are interested in the situation where the components 
have the same frequency but differ in polarization. 

By following the procedures outlined in Refs. [Tol l24j| . 
assuming only nearest neighbor interactions, and by us- 
ing the slowly varying envelope approximation, one may 
obtain the following model equations: 
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where L u u n = (C - Ak)u n + Vo(u n +i + u n -i) and 
L v v n = [C v n + Vo(v n +i + v n -i)]/si. Here f is the nor- 
malized propagation coordinate (y in the experiment). 
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The normalized envelopes u n and v n correspond to trans- 
verse electrically (TE) and magnetically (TM) polarized 
fields, respectively. The parameter (3 represents the non- 
linear coefficient, the normalized coupling constant is 
denoted by Vb, Afc is the normalized difference of TE 
and TM wave numbers, whereas Co can be regarded as 
a normalized propagation constant. By defining bire- 
fringence An = n x — n z and average refractive index 
= (n x -\-n z )/2, we write the function sj « 1+jAn/no- 
Finally, A = r xxz /r zzz and B = r xzx /r zzz , where r ijk 
denote the components of the Pockels tensor [lOj. One 
notice that in general, alike DVS in Kerr media 
17 , [1^ . in our situation there exists no possibil- 



ity to separate cross-phase and four- wave mixing effects. 
A conserved quantity of this model is the total power, 
P = E n (K| 2 + SiM 2 ) =P u + s 1 P v . By using Q, it 
can be shown that dP/dt; — implies: 

-_ e - OR, rV Im K<) 
^ ~ 81 d£ - l ^ B 2^ {l + \ Un? + sAvn? y 

(2) 

This expression gives us information on the energy 
(power) exchange among the TE and TM components 
which clearly will depend on the total level of power 
in each waveguide. By considering a one-channel 
constant-amplitude propagation of the form uq(£) = 
u exp[i(X u ^ + (j) u )\ and v (f) = v exp[i(X v ^ + V )], 
where A^ and & correspond to the respective propaga- 
tion constants and initial phases, respectively, we ob- 
tain the following expression for the power transfer: 
dP u /d^ ~ sin(AA£ + A</>), where A A = X v — X u and 
Acj) = (j) v — 4> u - By assuming only a linear and local de- 
pendence of the propagation constants, from (pQ) we get: 
A n ~ Co — Ak,X v « Co/si, which results in AA « Afc. 
If the components are initially in phase (Acj) = 0), the 
power transfer will be initially towards the TE mode pro- 
vided that Afc > 0, and towards the TM mode otherwise 

To gain a theoretical background of the model (pQ), we 
use a Newton-Raphson method to find coupled localized 
stationary solutions of the form: u n (£) = ii n exp(zA M £) 
and v n (£) = v n exp(zA v £), with u ni v n G 5ft. For the 
sake of simplicity we assume A n = X v = A, which in 
turn disables the power exchange between DVS compo- 
nents [Im(v n ii*) = in (|2j)]. These solutions may be 
regarded as the final stage of mode profiles after the 
DVS is formed. The power dependence of the two com- 
ponents on propagation constant, for the chosen set of 
experimentally achievable parameters [25], is shown in 
Fig. [2(a) in a logarithmic scale. The region of exis- 
tence of localized modes is between the low-amplitude 
and high- amplitude limits for the upper band edge plane 
wave [22| of the composed system of equations (pp). In the 
present case, this region corresponds to ~ A G { — 13,2}. 
Power of the TE mode always exceeds that of TM po- 
larization and, interestingly, grows in a similar fashion 
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FIG. 1: (Color online), (a) TE (dashed line) and TM (con- 
tinuous line) powers as a function of the propagation con- 
stant. Vertical lines separate regions, (b) and (c) TE and 
TM profiles, respectively: Diamonds, squares, triangles and 
stars correspond to A = —5,-1,0, and 1.1, respectively. 
Cq = Afc = p = 10, Vb = 1. 



as the power of the on-site mode A in Ref. 21]. For 
any A, the TE mode is always a one- hump structure [see 
Fig. njb)] which spreads transversally in the region of 
saturation [2l|. We may separate the PA diagram in 
smaller regions depending on the shape of the TM mode. 
In region I [~ A G { — 13,-0.6}], the TM mode corre- 
sponds to a one-hump structure [diamonds and squares 
in Fig. [11(c)]. In region II [-AG {-0.6, 1}], the TM mode 
corresponds to a two-hump structure separated by only 
one site [triangles in Fig. [2(c)]. In the next regions, the 
TM mode increases its distance between the two humps 
in an odd number of waveguides [as an example, see stars 
in Fig. ffl(c)]. 

While the total power increases, local saturation takes 
place [22]. As Afc > 0, the TE mode is the one which 
gains power. Therefore the TE mode starts to increase 
its power locally together with the TM mode [region I in 
Fig. UK a), diamonds and squares in Fig.[]](b,c)]. However, 
above some critical level of power [region II in Fig. [2(a), 
triangles in Fig. [2(b,c)] the local power in the center site 
is too high and the only possibility for the TM mode 
to exist is by exploring the neighborhood looking for a 
more stable configuration. Then, the TE mode further 
increases its power but now, due to saturation, by in- 
creasing the amplitudes in the next sites [see stars in 
Fig. [2(b)]. Again, the TM mode finds a new configu- 
ration which is initially stable, but now the separation 
between peaks consists of three sites [stars in Fig. [2(c)]. 
If we continue increasing the power we observe that the 
TE mode preserves its one- hump structure, by increas- 
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FIG. 2: Scheme of the experimental setup. 

ing its width, while the TM mode has a two-hump struc- 
ture where the separation between peaks continuously 
increases. Therefore, the DVS is mostly TE polarized, 
except at tails which have a dominating TM polariza- 
tion. The linear stability analysis of solutions coincides 
with the Vakhitov-Kolokolov criterion [5[ : modes are sta- 
ble for dP/dX > 0, and unstable otherwise. This implies 
that in region I solutions are always stable and, in the 
next regions, there exist both stable and unstable sub- 
regions [see Fig. QJa)]. 

To verify our theoretical predictions we use the ex- 
perimental setup sketched in Fig. El A cw laser with 
wavelength 532 nm is split into two orthogonally polar- 
ized (TE and TM) mutually coherent waves with the help 
of a polarizing beam splitter PBS. Optionally, to allow 
for mutually incoherent interaction of the two compo- 
nents, a TM polarized wave can be provided by a second 
laser of the same wavelength. Input power is adjusted 
with a combination of half wave plate A/2 and polarizer 
P. The two input beams are used to excite narrow single- 
channel TE and TM polarized modes of the WA by using 
a 40 x microscope lens ML. This nonlinear WA is fabri- 
cated in x-cut lithium niobate doped with copper. The 
length of our sample along the propagation y-direction 
is 11mm and the array consists of 250 parallel titanium 
in-diffused waveguide channels that are 4 jam wide with a 
separation of 4.4 fim (grating period A = 8.4 /xm) jifj. A 
second microscope lens ML images the intensity on the 
output face onto a CCD camera, where an additional po- 
larizer P allows for independent observation of both TE 
and TM components of the DVS. 

The nonlinear dynamics of TE only, TM only, and 
both TE and TM modes (mutually coherent from the 
same light source) is presented in Fig. 3. Here we make 
use of the fact that in photorefractives the nonlinearity 
grows exponentially in time, (3(t) = 0(1 — exp[— t/r]), 
where r is the dielectric response time [27]. Initially, af- 
ter switching on the light (t = 0), discrete diffraction is 
monitored for each situation. Although one may observe 
initial focusing of the TE mode within the first minutes in 
Fig. 03(a) (an even weaker effect is observed for TM), both 
modes alone are incapable to form a localized structure 
[Fig. 09(a,b)]. However, when both input polarizations 
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FIG. 3: (Color online). Experimentally observed DVS for mu- 
tually coherent input beams. Temporal nonlinear evolution 
of (a) TE component alone (P u = 150 /xW); (b) TM compo- 
nent alone (P v = 300 /iW); and (c) both components together 
{P u = 150 /iW, P. = 300 /iW) . 

are present [Fig. 09(c)] a five-channel wide DVS is formed 
after t = 30 min and remains stable for longer times t. 
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FIG. 4: (Color online). Mode analysis of stationary DVS 
solitons. (a) Stationary output (total intensity TE + TM) 
of DVS for different power ratios P u /P v ] (b) stationary TE 
and TM (amplified 40 x) polarized components for P u /P v — 
1; and (c) stationary TE and TM polarized components for 
Pu/Pv = 1 when input TE beam is blocked after formation 
of the DVS. 

Stationary images of DVS collected from the output 
facet of the sample for a fixed value of TE power and 
different values of TM power are presented in Fig. Ufa). 
As can be seen, the shape of the DVS slightly changes 
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for different power ratios P u /P v . TE and TM polar- 
ized components for a power ratio P u /P v = 1 (steady 
state) are shown in Fig. |U[b). As predicted, a dominat- 
ing single-hump TE polarized component and a weaker 
double-humped TM component are observed. The role 
of the TM input polarization can be further analyzed by 
blocking the TE input after stable formation of the DVS 
in Fig. [U(c). Obviously the TM polarized input light 
transfers most of its energy to the TE component, form- 
ing a single-hump solution, while the remaining power 
is trapped in form of a two-hump solution. This en- 
ergy transfer, from ordinary to extraordinary polariza- 
tion (TM — ► TE), is due to a specific anisotropic non- 
linearity in LiNbOs [in model (pQ), this corresponds to 
consider Ak > 0] . The mechanism of coupling of orthogo- 
nally polarized modes is explained by writing holographic 
gratings due to photovoltaic currents. Light is anisotrop- 
ically diffracted from these shifted gratings with polar- 
ization conversion, which leads to an energy exchange 
among the modes [28|. 
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FIG. 5: (Color online). Experimentally observed DVS for mu- 
tually incoherent input beams, (a) Stationary output of DVS 
[total intensity TE + TM] for P u /P v = 1.5; (b) TE polar- 
ized component; and (c) TM polarized components (amplified 
8x). 

Energy coupling of orthogonally polarized waves can 
be prevented by using mutually incoherent input beams 
(B = 0). Experimentally this is realized by a second laser 
of the same wavelength (see Fig. [2]), which now provides 
the TM polarized input beam, and corresponding results 
for the steady-state DVS formation are shown in Fig. [5l 
Again a two-hump structure is observed for the TM com- 
ponent, which now guides a significant part of the total 
power of the soliton. 

In conclusion, we suggested a rather general theoretical 
model to describe saturable discrete vector solitons hav- 
ing orthogonally polarized components. Power transfer 
and coupling between TE and TM components is inves- 
tigated as well as the corresponding localized stationary 
solutions. We discovered that these composite solitons 
might have different width and shape depending on the 
region of parameters. The dominating TE mode is single- 
humped while the weaker TM mode may exhibit both 
one- and two- hump structures. We confirm our findings 
experimentally by using either coherent or mutually in- 
coherent excitations, where the latter is used to suppress 



energy coupling in formation of discrete vector solitons. 
Our experimental conditions match the region II of sta- 
tionary solutions, a region with an intermediate level of 
power and highly localized solutions. This is because a 
one-channel input excites a strongly localized region of 
the array with high local intensity. The results obtained 
here could be useful in the codification of signals, filtered 
by polarization, in future all-optical communication net- 
works. 
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